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The unity of instantaneous spectral 
moments and physical moments 

By Jonathan M. Lilly 
NorthWest Research Associates, PO Box 3027, Bellevue, WA 98009, USA 

A modulated oscillation in two or three dimensions can be represented as the trajec- 
tory traced out in space by a particle orbiting an ellipse, the properties of which vary 
as a function of time. Generalizing ideas from signal analysis, the signal variability 
can be described in terms of kinematic quantities, the instantaneous moments, that 
formalize our intuitive notions of time- varying frequency and amplitude. On the 
other hand, if we observed an ellipse evolving in space we would seek to describe it 
in terms of its physical moments, such as angular momentum, moment of inertia, 
etc. The main result of this paper is to show that the two sets of moments are 
identical. Most significantly, an essential physical quantity — the circulation — is the 
same as the product of the two most important kinematic quantities, the instanta- 
neous frequency and the squared instantaneous amplitude. In addition to providing 
a rich set of geometric tools for the analysis of nonstationary oscillations in two 
or three dimensions, this result also has implications for the practical problem of 
inferring physical ellipse parameters from the trajectory of a single particle on the 
ellipse periphery, as is frequently encountered in the study of vortex motions. 

Keywords: Instantaneous frequency, instantaneous bandwidth, nonstationary 
signal analysis, trivariate signal, angular momentum, circulation. 

1. Introduction 

A powerful means emerging from the signal processing community for the analysis 
of nonstationary sign als is the what may be termed the method of instan taneous mo- 
ments, introduced by Gabor J 1946h and ex t ended by many author s (e.g. IVille L 1948 ; 
Bedrosiarl Il963t IVakman and Vainshteinl. 1977 ; Cohen and Led. 1 19891: iBoashasb , 
1992t ICohenl . Il995t IVakmanl Il996l IPicinbonol Il997t lLoughlin and Davidson! . l2000h . 
The instantaneous moments are functions of time, derived from the signal itself, 
that offer time-varying generalizations of our standard intuitive concepts — rooted 
in Fourier analysis, hence constant in time — of signal amplitude, frequency, and 
bandwidth. Since the instantaneous moments time-average to the corresponding 
Fourier-domain moments of the signal's spectrum, these quantities provide a direct 
and compelling link between time variabili ty and frequency-domain structure. This 
method had its genesis when lGaborl (Il94fih brought ideas from quantum mechanics 
to bear on outstanding questions in signal processing questions, formalizing the 
analysis of the information content of signals in a way that had not previously been 
possible. 

In recent years, as the intuitive conten t of the instant a neous moments has be - 
come more cl ear — t hanks to work such as Vakman ( 1996h . Cohen and Lee (1989), 
Cohenl (1995), and Louehlin and Davidson ~ (|2000l) — the virtue of this method for 
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the analysis of time-varying sig nals in observational data has b ecome increasing! 



appa r ent. In both geophysical (ITaner et all Il979t iRene et all 1986t iDiallo et al 



20061: iLillvL 1201 lh and oceanographic ([Lilly and GascardL l2006t iLillv et all 1201 ll 



applications, ideas based upon the method of instantaneous moments are proving to 
be remarkably informative for studying oscillatory processes that evolve with time. 
Thus ideas from physics brought to, and evolved in, the signal processing com- 
munity are now moving back into the physical sciences, providing a much- needed 
mathematical foundation for the analysis of nonstationary signals. 

A full appreciation of the power and potential of the method of instantaneous 
moments to the analysis of physical signals requires another generalization. This 
method has been developed for a univariate signal, i.e. a scalar-valued function of 
time x(t), whereas in the physical sciences, it is more frequently the case that the 
signal of interest is a vector x(t) in two or three dimensions. This is so because oscil- 
lations are typically generated by waves or wavelike processe s in phy sical space. A 
step in this direction was recently taken by ILillv and Olhede ( 2010rJl , who derived 
scalar- valued functions of time, called the joint instantaneous moments, that play 
the same roles for a multivariate signal x(£) that the standard instantaneous mo- 
ments play for a univariate signal x(t). The joint instantaneous moments have been 
shown to have informative expressions in terms of the geometry of the time- v arying 
ellips e traced out by the signal in two (ILillv and Olhede, 2010b) or three (jLillvL 



20111 ) dimensions. Recovery of multivariate oscillations from potentially noisy time 



series, together with estimation of the joint i nstantaneous moments, can then be ac- 
complished by a multivar iate generaliza t ion (|Lillv and Olhedd . [2012h of the wavelet 



ridge analysis method of lDelprat et al.l (|1992h 



As with the instanta neous moments of univaria te signal, the joint instantaneous 
moments introduced bv ILillv and Olhedd (2010b) constitute a scalar- valued func- 
tion of time at each order, corresponding to the orders of frequency-domain mo- 
ments of the signal's spectrum. Yet since the spectrum of a vector- valued signal is 
matrix- valued, we infer that the complete generalizations of instantaneous moments 
to such a signal should also be matrix-valued functions of time, and therefore the 
joint instantaneous moments must therefore only contain a fraction of the avail- 
able information. Such matrix-valued instantaneous moments have not yet been 
defined, and consequently the meaning of the additional information they contain 
has not yet been explored. In addition to achieving a more complete understanding 
of signal variability, there is a second compelling reason to examine in detail the 
instantaneous moment matrices: to close the apparent gap between the instanta- 
neous spectral moments on the one hand, and familiar physical moments on the 
other. The instantaneous moments could be aptly described as kinematic quanti- 
ties, in that they make no reference to physical properties. Yet if we have a particle 
moving in space we would immediately calculate its angular momentum as a func- 
tion of time, which is a physical moment. The link between physical moments and 
kinematic moments, if any, is not at all apparent. It is essential to establish this 
relationship if the method of instantaneous moments is to be used to illuminate the 
analysis of physical phenomena. 

The purpose of this paper is to fully generalize the instantaneous spectral mo- 
ments to the trivariate case, and to establish the connection between these quan- 
tities and standard physical quantities, building on work on the s calar- valued joint 
instantaneous moments for the bivariate and trivariate cases bv ILillv and Olhedd 
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(|2010bh and iLillvl (| 2 1 lh respectively. It is seen that the connection is surprisingly 
deep: the matrix- valued instantaneous spectral moments, and the physical moments 
of an elliptical ring of particles that is implicitly assigned to the time series by the 
foundational step in the analysis method, are identical. Most significantly, arguably 
the most important physical moment of the conjectured system — its circulation or 
average vorticity — is found to be identical to the product of the two most impor- 
tant kinematic moments, the instantaneous frequency and squared instantaneous 
amplitude. This has significant implications for the recovery of vortex properties 
from particle trajectories , a practical proble m that is encountered in oceanography 
(jLillv and Gascardl l2006t iLillv et all [201 lh and in fluid dynamics more generally. 

Fundamentals of the method of instantaneous moments for multivariate signals, 
together with necessary preliminary material from the signal processing literature, 
are presented in ^ The main result — the unity of moments — is stated without 
proof in [21 an d its implications are discussed. The proof of the unity of moments is 
given in []4l In SJSJ the interpretations of the instantaneous moments are investigated 
in more detail. For reasons of space, this paper focuses on theoretical as opposed to 
practical issues. The connection to practical applications is briefly discussed along 
with the conclusions in ^51 



2. Fundamentals 

The variability of a trivariate signal can be usefully described in terms of a parti- 
cle orbiting a time- varying ellipse, building on a powerful method from the signal 
processing literature for linking time-domain variability to frequency-domain struc- 
ture. In the trivariate case, frequency-domain structure is captured by the 3x3 
spectral matrix. Recent work has established the relationship between the ellipse 
geometry and the trace of this matrix, but the connection to the remaining terms 
has not yet been investigated. 



(a) An elliptical ring of particiles 

A real- valued signal x(t) = \x(t) y(t) is observed in three dimensions, 

which for clarity is imagined to be the trajectory traced out by a hypothetical 
particle in physical space. The signal is regarded as deterministic and zero-mean, 
and is assumed to have been observed over its entire duration and to be of finite 
energy, i.e. \\x(t)\\ 2 dt < oo. Here ||z|| 2 = z H z is the vector norm of some 
possibly complex-valued vector z, with the superscript "H" denoting the Hcrmitian 
transpose. There is reason to believe that the trivariate signal x(i) is oscillatory in 
nature, but with properties that may change in time. Our goal is to analyse x(f) in 
such a way as to illuminate its structure. 

The kinematic model of a time- varying ellipse will form the starting point. The 
evolving position of an elliptical ring of particles, constrained to remain coplanar 
and elliptical in shape but otherwise free to evolve in three dimensions about its 
centre of mass, may be written as 



x(i )¥ >) = Q x (i) 



a(t) ' 
b(t) 




cos 4>(t) — sincj>(t) 
sin 4>(t) cos (f>(t) 



COS(/9 

sin ip 



(2.1) 
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where ip is a particle label with — ir < ip < 7r. Here a(t) and b(t) are the time- 
varying ellipse semi-major and semi-minor axis lengths with a(t) > > 0, and 
the orbital phase <p(t) specifies a time-varying shift of particle locations around the 
ellipse periphery. Q x (t) is a 3 x 3 rotation matrix written as 

Q x (i) = 3 3 (a(t))J 1 (/3(t))3 3 (e(t)) (2.2) 

and represents a time-dependent rotation in the so-called "Z-X-Z" form, where 



Ji(/3) 



10 

cos (3 — sin f3 
sin j3 cos (3 



J3(«) 



cos a — sin a 
sin a cos a 
1 



(2.3) 



are rotation matrices about the x and z axes, respectively. A sketch of the ellipse 
geometry is shown in figureQ] In this rotation, 8(t) sets the orientation of the major 
axis of the ellipse in a two-dimensional plane, while f3(t) and a(i) are the zenith 
and azimuth angles of the normal to the plane containing the ellipse. From right to 
left, the kinematic model (|2.ip states that a particle labeled by phase ip is shifted 
4>{t) radians counterclockwise around the unit circle in the x-y plane, which is then 
deformed along the x-axis and y-axis to generate an ellipse, and which is finally 
subjected to a rotation in three dimensions. As t varies with fixed ip, x(t, ip) traces 
out the trajectory in three-dimensional space of the particle labeled by ip. 

This kinematic model has considerable dynamical relevance. Setting a = f3 = 
leads to a two-dimensional version of (|2.ip that emerges frequently in fluid dynam- 
ics, encompassing the oscilla ting ell i ptical vort ex solution t o the rotating shallow 
water equations e xamin e d by | Younel l)l98fih and lHolml fll99lh as well as the archety- 
pal Kirchhoff and iKidal (|198ll ) family of two-dimensiona l vortex solu tions. In fact, 
(|2.1I) with a = f3 = is identical with equation (2.6) of iHolml (| 1 9 9 lh after substi- 
tuting his (4.2). In that work, a set of evolution equations for the ellipse parameters 
is derived in which periodic motion arises as a consequence of physical conservation 
laws. For this and other vortex solutions, the kinematic model (|2.ip would be seen 
as generating a family of trajectories on the left-hand side based on prescribed or 
physically determined ellipse parameters on the right-hand side. The primary con- 
cern here will be the opposite, namely, given a signal we will seek to assign or infer 
a set of ellipse parameters in some meaningful way, thus moving from the left-hand 
side of (|2.1[) to the right-hand side. 



(b) The canonical ellipse parameters 

An infinite family of ellipse parameters can generate the same physical signal. 
The method for assigning a unique set of physically meaningful time- varying ellipse 
parameters to a trivariate signal x(t) beings with the construction of the analytic 
signal, the starting point for the metho d of instantaneous momen ts. H ere w e follow 
the recent mu l tivari ate d evelopment in iLi lv and Olhedel (l2010bh and lLillvl (l2f)llh : 
see iBoashashl (|l992t ) and Picinbonol ( 1997 ) and references therein for background. 
With Hx.(t) = i du being the Hilbert transform of x(£), where "f" is the 

Cauchy principal value integral, the analytic part of x(£) is given by 



x+(t) = 2Ax(t) = x{t) + iHx(t) = e i<#,(t) Q x (i) 



a(t) ■ 
-ib(t) 




(2.4) 
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Ellipse Schematic in 3D 




X-Axis Y-Axis 



Figure 1. Schematic of an elliptical ring of particles in three dimensions. The ellipse has semi-major 
axis a = 3, semi-minor axis b = 3/2, precession angle 9 = it/3, zenith angle j3 = 7r/4, and azimuth 
angle a = n/G. The phase increases from an initial value at <j> = 5n/G, tracing out the heavy back 
curve through one full cycle, during which time all other ellipse parameters are constant. The 
plane of the ellipse is indicated by the dotted lines, with the original x- and y-axes marked by 
thin dashed lines. A heavy dashed black line shows position of the "particle" at the initial time, 
while a heavy gray dashed line marks the ellipse semi-major axis. A heavy solid black line, with 
a black circle at its end, is the normal vector to the plane of ellipse; the projection of this vector 
onto the x— y plane is shown with a heavy solid gray line. 



where A is called the analytic operator. This defines the six ellipse parameters on 
the r i ght-ha nd side in terms of the three complex- valued quantities in x_|_(i); see 
iLillvl (|201lh for explicit expressions. Comparison with the kinematic model shows 
that (|2.4[) corresponds to setting <p = in (|2.1[) . thus x(i) = x(t, 0). Thus when 
a signal x(£) is observed, this is regarded as representing the position of one of a 
continuum of particles in an imaginary elliptical ring, the properties of which are 
assigned to the signal by the action of taking the analytic part. We refer to this 
model for the evolution of x(t) as the modulated ellipse representation. 

The ellipse parameters assigned to x(t) by t he analytic signal method are called 
the ca nonical ellipse parameters. It is argued in iLillv and O lhcde (2010b) and|Lillv 
([20111 ). for the bivariate and trivariate cases respectively, that the ellipse parame- 
ters so determined represent the natural or intrinsic choice from the infinite family 
of possibilities. The modulated ellipse representation explicitly captures joint vari- 
ability by considering the components of x(i) to be aspects of a single object. This 
can be seen as a natural time- varying generalization of the representation of the 
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Fourier spectrum of a trivariate signal i n terms of ellipse geom etry as a function of 
frequency, as is common in optics (e.g. Born and Woll . 1970l) . oceano graphy (e.g 



Caiman . 19781) . and seismology (e.g. lPark et al. . 19871) . for example. See iLillvl (|201lh 



for further discussion of the ubiquity of the ellipse perspective in different fields. 

(c) The analytic spectral matrix 

The aggregate second-order structure of x + (i) is contained within its spectral 
matrix, or its Fourier transform, the autocovariance function. If X(w) is the Fourier 
transform of x(i), with x(t) — X(w)e lwt do;, then we may define the one- 

sided deterministic spectral matrix as 

s x h = x+hx^h (2.5) 

where X+(w) = 2[/(o>)X(cj) is the Fourier transform of x+(£) with U{oj) being the 
unit step function. Frequency-domain structure is then conveniently analysed in 
terms of the global moments of the spectral matrix 

S n = — — f uj n S x (cu)duj (2.6) 
27T Jo 

which are themselves matrix- valued. The Fourier transform of the spectral matrix is 
the deterministic autocovariance function associated with the analytic signal x + (t) 

/OO -i />oo 

x+(i + r)x£(t)dt = — / S x (wy^duj. (2.7) 

An expansion of the autocorrelation matrix in terms of the time lag r then gives 

Rx(r) = So + irSi - \r 2 S 2 + ■■■ (2.8) 



in which the frequency-domain moments defined in (J2T6J) appear as the coefficients. 
Thus the global moment matrices S„ can equivalently be seen as describing as- 
pects of the aggregate frequency-domain structure of S x (w), or else aspects of the 
aggregate time-domain structure captured in R x (t) via powers of the time lag r. 
For future reference, we introduce S x (w) = trS x (w) = ||X(w)|| 2 as the trace of the 
spectral matrix. This scalar-valued function of frequency, called the joint analytic 
spectrum, is the sum of the spectra of the analytic parts of all signal c ompon ents. 
The Joint a nalytic spectrum S x (lu) was examined by iLillv and Olhedc (l2010bh and 
Lillvl (|201lh . and this work will make an important reference point in what follows. 



The use of the one-sided Fourier spectrum is standard in many fields, and it is 
important to emphasize that this involves the implicit use of the analytic signal. 
Two crucial aspects are however not addressed by analysis of S x (cj) or R x (r). The 
first is that there is not yet a facility for investigating time-varying structure. This 
is significant because many interesting signals exhibit properties that evolve with 
time, yet can have spectra which are indistinguishable from other, completely dif- 
ferent signals. The second is that is not obvious how to assign an interpretation to 
the moments S„. One could discuss coherences, for example, but such an approach 
describes relationships between pairs of channels, rather than addressing the sig- 
nal as a unified object with structure extending across all channels. This suggests 
connecting the matrix- valued Fourier-domain moments S„ to the modulated ellipse 
model for time-varying structure, which is to be accomplished herein. 
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(d) The joint instantaneous moments 

The analytic signal is the basis for the creation of key time-varying quanti- 
ties, called the instantaneous moments, that allow the moments of the spectrum 
of a univariate s ignal to be decomposed across time (e.g. iGabor , Il946t IVirleUl948t 
Cohen and LeeL Il989h . In particular, the first three of these quantities — the in- 
stantaneous amplitude, frequency, and bandwidth — provide a pow erful link between 
time va riability and frequency-domain structure. Here we follow iLillv and Olhedc 
(|2010bl ) in extending these quantities to describe the scalar- valued joint analytic 
spectrum S x (u) = trS x (w) of a multivariate signal. Introduce 



4(*) = -l|x + (t)|| 2 



Im{x«(t)x;(t)} 



l|x+(t)| 



x+(t)|| 2 

2 



«x(<) 



(2.9) 
(2.10) 

(2.11) 



called the joint instantaneous amplitude, frequency, and bandwidth, respectively, of 
the multivariate signal x(i). These are found to satisfy the integral constraints 



r°° i i f°° i 

J 4(t)dt=-—j S x Hd,E-£ 

1 [°° 1 f°° 

., / w x (£)||x+(i)|| 2 d£ = — — / ujS x (uj)duj = uj 

£x J-oc 27r t x J 



(2.12) 
(2.13) 



together with 



1 f°° r i 
-J^ {vl(t) + -EJ X ] 2 } ||x + (t)|| 2 dt 

1 f°° 

^7T^/ k^xfSxMdwE?^ (2.14) 

and thus K x (t), u> x (t), and u x (i) respectively give the essential instantaneous or 
time-varying contributions to the zeroth-order, first-order, and second-order mo- 
ments of the joint analytic spectrum S x (uj). In the above, £ x /2 is the total energy 
of the signal x(t), uj x is the mean frequency of S x (oj), and a x is the standard de- 
viation about the mean frequency. The name "instantaneous bandwidth" for v x (t) 
arises from the fa ct that W x is a measure o f the signal's global Fourier bandwidth. 

It is shown in ILillv and Qlhe"del ( 2010bl ) that the joint instantaneous moments 
defined by (|2.9[) - (|2.1ip . reduce, for a univariate signal, to the standard definitions 
of the instantaneous amplitude, frequency, and bandwidth, respectively. Therefore, 
as argued therein, the joint instantaneous moments generalize the univariate in- 
stantaneous moments to accommodate the moments of the joint analytic spectrum 
S x (uj) of a multivariate signal. The integral relations (I2.12p — (12.141) a r e preci sely 



those which occur in the univariate case, see §11 B of ILillv and Olhedd (l2010bh . 

As discussed in more detail later, the joint instantaneous moments n x (t), uj x (t), 
and v x (t) can be shown to describe the variability of the oscillatory signal x(£) in 
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terms of the geometr y of the modulated ellipse repr esenta t ion. T his interpretation 
has been pursued by iLillv and Olhedd (|2010bh and iLillvl ([20111 ) for the bivariate 
and trivariate cases, respectively. However, in the trivariate case, the Hermitian 
instantaneous moment matrices S n (£) each contain nine independent quantities, 
three real- valued quantities along the diagonal and three complex- valued quantities 
in the off-diagonal terms. The joint instantaneous moments, being the traces of 
these matrices, describes only one of the nine components. In the following we will 
further generalize the instantaneous moments to matrix-valued quantities, and then 
obtain simple physical interpretations for the remaining eight terms. 



3. Unity of kinematic and physical moments 

Matrix- valued instantaneous spectral moments — kinematic quantities which gener- 
alize signal processing concepts of instantaneous amplitude, frequency, and band- 
width to a multivariate signal — are introduced. These are found to be identical to 
the hierarchy of physical moments of the evolving elliptical ring of particles which 
uniquely describes the observed signal, with several notable implications. 



(a) Matrix-valued instantaneous moments 

The analytic signal method described in the previous section can be used as the 
basis for achieving a time-varying decomposition of the spectral matrix S x (a;) of 
the sign al x(t), building on the joint instantaneous moments of Lilly and Olhedei 



(|2010bf ) that decompose the matrix trace. Define 

S (t) = *+(*)*+(*) (3-1) 
Si(t) = hcrm{-ix' + (i)x^(£)} (3.2) 

S 2 (t) = (3.3) 

as the zeroth-order, first-order, and second-order instantaneous moment matrices, 
respectively, where hermF = (F + F H )/2 is the Hermitian part of the matrix F. 
Observe that these satisfy the integral constraint 

-i /"OO />oo 

S„ = — / u n S x (uj)du= S n (t)dt (3.4) 
2tt Jo J_oo 

and consequently the S n (t) give the instantaneous or local contributions that inte- 
grate to the corresponding global Fourier-domain moment matrices S n . The S n (t) 
are not uniquely defined by the integral condition (|3.4|l. However it i s sensible, and 
conventional in the analysis of univariate sig nals fe.g. ICohenl . Il995h . to constrain 



the instantaneous moments to have Hermitian symmetry, like the global moments 
to which they integrate; this leads to (|3.1|) - (|3.3|) as the natural definitions. 

The quantities So(i), Si(i), and S2(t) provide normalized matrix-valued gener- 
alizations, appropriate for a multivariate signal, of standard instantaneous moment 
quantities of a univariate signal — the instantaneous amplitude, frequency, and sec- 
ond central moment, respectively. As these are all Hermitian matrices, each contains 
nine real- valued functions of time. Taking the trace of (|3.4I) . we find 



/£(*) = |tr{S (t)}, = V ^ = J^W-^- ^ 
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and the traces of the matrix-valued i nstantaneous moments re cover the scalar- 
valued joint instantaneous moments of iLillv and Olhedl (|2010rA . The remaining 



terms in the S n (t) capture the time variability that accounts for the full second- 
order frequency-domain structure implicit in the spectral matrix S x (w). Note that 
in defining the instantaneous moment matrices, it proves convenient to eschew the 
convention of normalizing by the signal power ||x + (t) || 2 , as is done in the univariate 
case and for the joint instantaneous moments. According to these definitions, the 
second instantaneous moment matrix S2(i) of x(i) is also the zeroth instantaneous 
moment matrix of a different signal, the time derivative x'(i) of the original signal. 

(b) Physical moments of an elliptical ring of particles 

Using the analytic signal, an observed trivariate signal x(t) can be described as 
if it were due to the trajectory traced out by a particle orbiting the periphery of 
a time-varying ellipse. The time-varying position of any other particle within this 
hypothetical ellipse may then be found by a suitable choice of phase shift. That is, 
we define a new signal x(i, ip) as 

Re{e iip x+(t)} = x(t,ip) (3.6) 

and observe that this recovers the form of the kinematic model (|2.1j) . By varying 
the phase shift 95, the trajectory of the observed particle x(t) — corresponding to 
x(t, 0) — can be transformed into the trajectory x(i, ip) of any other particle in the 
ellipse. The modulated ellipse representation of a signal may therefore be viewed 
as specifying the time-varying position of an entire elliptical ring of particles, the 
canonical ellipse. This suggests computing the physical moments of the canonical 
ellipse, in order to understand how these are related to the kinematic moments. 

Some physical moments of a ring of particles lying in a plane are defined in 
Table [T] The ring is taken to have mass density p(x) = dM/d£, with a total mass 
M = <f> c /o(x) d£ equal to unity, and with its centre of mass at the origin. The mo- 
ments are sorted in terms of their order — that is, the number of time derivatives 
involved — as well as whether they are scalar-valued, vector- valued, or tensor- valued. 
There are two kinds of integrals in Table [TJ integrals in which the variable of in- 
tegration is the scalar-valued differential arc length d£, occurring in the moments 
that lie along the main diagonal and counter diagonal of the tabic, and directed 
line integrals with a vector- valued variable of integration dx. Moments of the latter 
type are defined to include a weighting by the linear mass density p(x) , while those 
of the former type are defined to include a factor of ^ but no mass weighting. It 
will emerge later that with the assumption of a unit-mass ring made here, both of 
these choices are equivalent; the distinction in Table Q] is made in order that the 
moments take on familiar forms. 

In Tableland henceforth, I is the 3x3 identity matrix. Also we have introduced 
a notation to map a vector f = \f x f y f z ] into a skew-symmetric matrix, 

" /, ■ 

erf = f z -f x . (3.7) 

-fy fx _ 

The notation "cr f " is a reminder that we may write the cross product f x g as the 
matrix multiplication (erf) g, and thus erf may be referred to as the cross-product 
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Table 1. Some moments of a ring of particles 





Scalar-valued 


Vector-valued 




Matrix-valued 











= 27fc X 


x dx 


= 4b (!| 


x| 


2 -xx T ) pdl 


1 


S - 2n §C U dX 


m x (t) sfcxx 


upd^ 


^T x (t) = herm 


i 

L>5T 


§ c (cr u) T cr dx 


2 


«i(t)=/ ||u|| 2 pdi 


n„(t) = £ f u 


x du 


2Tu(t) = / (I| 


u 


2 -uu T )pd^ 



(Moments of a time- varying closed curve of particles C(t) having unit mass, with p = p(x) 
being the linear mass density along the curve. The particle position is x and the velocity 
is u = u(x). Directed line integrals with a vector- valued variable of integration are taken 
by convention in the right-hand sense. The moment order is shown at the left.) 



matrix associated with f . The cross-product matrix is in use in various communities 



e.g. Nour-Omid and Rankinl . ll991tlLiu and Trenkleii 12008: Baks alarv and Trenkler 



20111 ). although a common notation does not seem to have yet emerged. 

The categorization in Table [T] includes many familiar moments as well as some 
less familiar ones. A measure of the typical distance from the ring to the origin 
is given by K x (t). The circulation r x (i) and kinetic energy | K„(f) are familiar 
physical properties, as are the spatially averaged angular momentum vector m x (t) 
and the moment of in ertia tensor X x (t). The vector n x (£), termed the normal vector 



following iLillv ( 2011 ). is more rarely encountered, although this integral appears in 



the c alculation of the magnetic moment of a current loop (e.g. Arfken and Weberl . 



19951 pp 64-65). This vector points in the direction normal to the ring of particles, 
and as discussed in the preceding reference, from a variant of Stokes' theorem its 
magnitude is found to be 1/ir times the area enclosed by any curve C. At second 
order, n u (t) plays a similar role in describing the direction that is normal to the 
velocity, while X u (t) is the kinetic energy tensor defined in such a way as to mimic 
the form of the moment of inertia tensor I x (t). The most unusual quantity is the 
first-order tensor J7" x (£) ; this is a matrix- valued mixed product between the velocity 
and the position, complementing the scalar- valued mixed product, the circulation 
r x (i), and the vector-valued mixed product, the angular momentum m x (i). 

The main result of this work is to show that the first three instantaneous mo- 
ments matrices, defined in p.l[) - fl3.3[l in order to identify the instantaneous contri- 
butions to the Fourier-domain moments of the spectral matrix, take the forms 

iso(t) - 4(t)I-iicrn x (i)-X x (i) (3.8) 

isr(t) = E^i + i l crmx ( t) _j- x(t) (3. 9 ) 

is 2 (i) = /£(t)I-iicm u (t)-Z u (t) (3.10) 

and consequently can be interpreted as containing the physical moments of the 
canonical ellipse. As will be shown later, each of these complex- valued Hermitian 
matrices contains, from left to right, a scalar-valued portion, a vector-valued por- 
tion, and a tensor-valued portion, respectively. The instantaneous spectral mo- 
ments, and the physical moments of the canonical ellipse, are therefore identical. It 
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is not at all obvious that this should be the case. This result is termed the unity of 
kinematic and physical moments. 

Some implications of this result are pursued in the remainder of this section. 
The proof of the unity of moments is postponed until Detailed interpretation of 
individual terms in the instantaneous moments will be addressed in Sj5] 

(c) Interpretation of instantaneous contributions to the spectral matrix 

The unity of moments provides a rich set of geometric tools for analyzing the 
time-varying structure of a bivariate or trivariate signal. The structure implicit 
in the matrix-valued instantaneous moments can now be described using a set of 
physical concepts, such as moment of inertia and angular momentum, with which 
we are already intimately familiar. This interpretation provides us with new insight 
into the spectral matrix itself. 

For example, the zeroth-order global moment So of the spectral matrix can be 
decomposed as 



where the averages on the first line are defined as the corresponding integrals over 
instantaneous quantities on the second line. The scalar, vector, and tensor portions 
of the zeroth-order instantaneous moment matrix are therefore seen as being re- 
spectively built up from instantaneous contributions from the squared instantaneous 
amplitude, the normal vector, and the moment of inertia tensor of the canonical 
ellipse. The variability of So(t) therefore captures, in its different structure com- 
ponents, aspects of the geometric variability of the canonical ellipse which imprint 
themselves directly onto the spectrum. This lends itself to a wide variety of analysis 
applications. One could, for example, temporally average So(i) and the other in- 
stantaneous moment matrices during different time periods in order to quantity the 
contributions to the spectral matrix due to distinct modes of observed behaviour. 

More generally we may see the scalar-valued, vector- valued, and tensor-valued 
portions of the instantaneous moment matrices can be seen as decomposing the 
corresponding portions of the spectral matrix itself, which can be written 



as will be shown in the next section. Observe that the quadrature spectrum — the 
imaginary part of S x (w) — is entirely wrapped up in the vector- valued term in (|3.12l) . 
Thus, from (|3.8[) - (|3.10l) . we see that the zeroth-order, first-order, and second-order 
instantaneous moments associated with the quadrature spectrum are the cross- 
product matrices of the normal vector of the canonical ellipse n x (i), the negative of 
the angular momentum vector m x (i), and the normal vector of the time derivative 
of the canonical ellipse n u (i), respectively. This gives a new way of interpreting 
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the information contained within the quadrature spectrum. Alternatively, the com- 
ponents of instantaneous moment matrices may be viewed as providing different 
contributions to the autocorrelation matrix R x (t) associated with various powers 
of time lag r, which we see from inserting (|3.8|) — (|3.10[) into (|2. 



(d) Recovery of physical ellipse parameters 

A second important implication concerns those signals for which the physical 
moments in Table [T] are in fact themselves the sought-after quantities. Imagine 
that there exists a physical ring of particles in space, evolving according to the 
kinematic ellipse model (12.11) for some choice of time-varying ellipse parameters. 
A broad class of physical systems could be described in this manner, including 
observations of oceanic vortices using freely drifti ng or "Lagrangian" instruments 
(|LiHv and Gascard . l 2006t iLillv and Olhedel . l2012h . In general, the trajectory of a 
particle in such a system will become associated, after we take the analytic part, 
with a set of canonical ellipse parameters that are not the same as the physical sys- 
tem that generated the signal. If we put an arbitrarily chosen set of time-varying 
ellipse parameters into the modulated ellipse representation (I2.4p . there is no rea- 
son to assume that this will generate an analytic signal, and yet we know that the 
canonical ellipse parameters will do so by definition. If on the other hand the phys- 
ical ellipse parameters are such that they do generate an analytic signal via (|2.4[) . 
then the analytic signal method will exactly recover them from observed signal. 

This means that in some circumstances, the trajectory of a single particle is 
sufficient to recover exactly the time-varying properties of the entire system, and 
these physical moments integrated across time determine the structure of the ob- 
served spectral matrix. Denote the physical ellipse parameters as a(t), b(t), etc., 
in order to emphasize that these are prescribed and not necessarily the canonical 
ellipse paramete rs of the generated signal. It follows immediately from the work of 
Bedrosianl (1963) for univariate signals — the so-called "Bedrosian's theorem" — that 
a sufficient condition for the physical ellipse parameters to be canonical is that all 
components of the 3x2 matrix 



Qx(t) 



a{t) J) ' 
b(t) 




:J 3 (S(t))Ji (/?(*)) J 3 (0(*)) 



a(t) J) ' 
b(t) 




(3.13) 



are supported entirely on lower (i.e., less positive) frequencies than is the phase 
function e 1 "^*' . This amounts to a statement that the ellipse geometry evolves more 
slowly than the circulation of the particles around the ellipse. 

Under these conditions, Bedrosian's theorem implies that the analytic operator 
"passes through" the matrix, giving 



x+(i) = 2A{ Q x (t) 



a(t) ' 
b(t) 




COS (j){t) 

sin <f)(t) 



Q*(f) 



a(t) ■ 
-ib(t) 




(3.14) 



and since this matches the form of (|2.4p . we see that the prescribed and canoni- 
cal ellipse parameters are identical. As it is reasonable to believe that a subset of 
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oceanic vortices may satisfy or approximately satisfy this condition of slow varia- 
tion, it appears that the link between the observed trajectory of a single Lagrangian 
particle, and physical structures in the vicinity of this observation point, is consid- 
erably more direct than has been previously recognized. Further investigation of 
this point, and comparison between the domains over which the ellipse parameters 
are canonical and the domains over which dynamical solutions exist, is called for 
but its outside the scope of this work. 



4. Proof of the unity of moments 

The proof of the unity of moments consists of three components: a Hermitian matrix 
expansion; a relationship between products of analytic signals and a phase averaging 
operator; and a reduction of integrals around the ellipse to a phase average. 



(a) A matrix expansion using complex vectors 

In what follows we will make use complex-valued three- vectors, which warrants 
some discussion as these objects are somewhat rarely encountered. However, their 
use dates back over a century to Gibbs and the beginn i ng of modern vector analysis 
(see "Note on bi vector analysis", p. 84 — 90 of iGibbsl. 119061). N ecessary ba sic facts 
are gi ven h ere; further details m a y be found in Gibbs " lfl906h . as well as Lindell 
fll988h and iBoulanger and Haved (Il993h . The cross-product of two real-valued or 



complex- valued vectors f = [f x f y f z ] and g = [. 



9x 9 V 9z\ 



is defined as 



f X g = {fy9z - fz9y) 1 - (fx9z - fz9x) j + (fx9y ~ fy9x) k 



(4.1) 



where i, j, and k are the unit vectors along the x, y, and z-axes, respectively. Observe 
that there is no conjugation in this definition. The cross product of a complex vector 
with itself vanishes, f x f = 0, but the cross product of a vector with its own complex 
conjugate has a nonzero purely imaginary part, f x f * = 2ilmf x Ref. 

The following identities involving the cross-product matrix of complex-valued 
vectors may be readily verified: 



(crf) H crg= (i H g)I- K f H 



— crlm {f x g*} 



Im {herm gf ff } 



(4.2) 
(4.3) 



where the cross-product matrix cr f of a real- valued or complex- valued vector f is 
de fined by (13.71). Other ide ntities involving the cross-product ma trix may be found 



Liu and Trenklerl (|2008) and Baksalarv and Trenklerl <|201lh . Combining (|4T2| 



and (|4.3I) yields an expansion for a Hermitian matrix, which we write without loss 
of generality as herm gf H , as 



(i) 



(ii) 



hermgf ff = Re |f ff g} I — i- cr Im {f x g*} — hermRe |(crf) H crg|. (4.4) 

This consists of a scalar- valued product (i), a vector- valued product (ii), and a 
matrix-valued product (iii) of f and g. The scalar-valued and matrix-valued portions 
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(i) and (iii) are real-valued, while the vector-valued portion (ii) is imaginary. Note 
that the traces of the three quantities on the right-hand side are 3 Re {f^g}, zero, 
and —2 Re {f H g}, respectively; thus the coefficient of I in (I4.4[) . which we refer to 
as the scalar-valued portion of the matrix, is also the matrix trace. 
Using (|4.4[) . we now expand the instantaneous spectral matrices as 

S (t) = ||x + (t)|| 2 I-iicrIm{x + (£) x <(*)} 

- hcrm Re j [cr x+ {t)] H cy x+ (t)\ (4.5) 

for the zeroth moment matrix and 

Si(t) =Im{x^(f)x' + (f)}l + iicrRe{x + (i) x x+{t)} 

- herm Im { [cr x+ (i)] H cr x' + (t)\ (4.6) 

for the first moment matrix, and so forth. In the following subsections we will sec 
how to relate these various terms to the physical moments. 



(b) Phase averaging 

Products of the components of the analytic signal x+(£) or its derivatives im- 
plicitly evaluate integrals around the periphery of the canonical ellipse, as will now 
be shown. Let f(t) be the product of some number, say M, components of x(i) or 
one of its derivatives of some order 

f{t) = x^\t)x^\t)---x ( t M M \t) (4.7) 

where the superscript "(n m )" indicates the n m th-order derivative with respect to 
time; here we assume derivatives of up to the required order exist and are finite. 
The phase average of /(t), defined as 



f{t) =i j v < 1 '(t^)< 2 '(t,V>)--^\T ) (t^)^ (4.8) 

produces a smoothed, but still time-varying, version of some property f(t) by av- 
eraging over all particles in the ellipse. The terms appearing in the integrand are 
components of the phase-shifted version of the signal x(t, ip) defined in (I3.6|) or one 
of its partial derivatives with respect to time, which we indicate with the super- 
scripts. Effectively, (|4.8p allows us to separates variation around the ellipse from 
ch anges in the ellipse geom etry; variation s on th is idea have been used previously 



bv lLillv and Gascardl (|2006l ) and ISchreierl (|2008l ). 



For a given function f(t) of the form (|4.7p . the phase average is constructed by 
algebraically calculating expressions for the M terms in f(t) using the modulated 
ellipse representation of the phase-shifted signal (|3.6I) . and then performing the 
indicated average in (I4.8[) . As an example of a quadratic product (M = 2), the 
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phase-averaged squared magnitude of the real-valued vector x(t) is given by 



h(t)P = KOI 2 +W0P +K0P =s= / x'ftrfxft^dy 

1 

— / [a 2 {t) cos 2 (4>(t) + (p)+ b 2 {t) sin 2 (<f>(t) + (p)] dip 



2tt 



a 2 (i) + 6 2 (i) 1„ ,. ||a 



~iM*)ir = <w (4-9) 



which is one-half the instantaneous power of the analytic signal vector. This does 
not depend on the orbital phase </>(£), while by contrast the unaveraged quantity 

||x(t)|| 2 = a 2 (t) cos 2 m)) + b 2 (t) sm 2 (<f>(tj) (4.10) 

varies rapidly as the particle orbits the ellipse, even for fixed ellipse geometry. 

In this way the analytic signal can be used to greatly simplify calculations 
involving phase averages of the real- valued signal. Let £i(t) and £2^) each be an 
element of x + (t) or one of its derivatives, so that Re £1 (t) and Re £ 2 (t) each represent 
one of the M terms in (|4.7I) . £1 (t) and £2 (t) are analytic because taking the analytic 
part commutes with differentiation. The phase averages of the products of real or 
imaginary parts of £i(i) and £2(2) are found to be 

Re £x (t ) Re £ 2 (t ) V = Im & (t ) Im £ 2 (t ) V = i Re {£j (t ) £ 2 * (i ) } (4.11) 

Im£ 1 (i)Re£ 2 (i) V = -Re£i(i) Im£ 2 (i) * = ~ Im {&(*)£(*)} (4.12) 

and therefore, significantly, these phase averages can be calculated by simply taking 
the real or imaginary part of the Hermitian product of the analytic functions. In 
the above example, direct calculation in (|4.9[) shows that ||x(i)|| 2 = ||x+(£) || 2 /2. 
However, this follows at once from (|4.11l) . because choosing £i(t) = £2(0 = x + (t) 
we have |^(i)| 2 = \x + (t)\ /2 and similarly for y(t) and z(t). 

To verify the quadratic phase averaging results (|4.11[) and (|4.12l) . write the 
functions £i(i) and £2(0 in the form 

£x(t) = [ci( t) + [dl (t)] , £ 2 (i) = e 1 ^) [ca(t) + id 2 (i)] (4.13) 

where c\ it) , d\ {€) , C2 (t) , and <i 2 (*) are all- real valued and do not depend explicitly 
on the phase <f>(i), although they may depend upon derivatives of 4>{t). To see (|4.11l) . 
note that the product of the real parts is given by 

Re £1 (t) Re £ 2 (t) - - [c x (i)c 2 (t) + d x (t)d 2 (t)] 

+ - [ci(i)c 2 (t) - di(t)da(t)] cos (2^(t)) 

- i [ Cl (t)d 2 (i) +c 2 (t)d 1 (i)] sin (20(i)) (4.14) 

where we have made use of a trigonometric identity. At the same time, we have 
Reld-MSM} = c i(t) c 2(t) + di(t)d 2 (f), and thus carrying out the phase average 
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in (|4.14[) obtains (|4.11l) . Similarly Im,£i(i) Im^it) is found to contain a constant 
term, equal to the constant term in (I4.14[) . plus oscillatory terms. We then find 

Im^Re^r = -Re6W Im6(t)" = \ \c%{t)d x {t) - Cl (t)d 2 (t)] 

= 5 Im {&(*)£(*)} (4.15) 

for the phase average of the product of real and imaginary parts, verifying (I4.12[) . 
In |Appendix A| a related phase averaging result for a quartic (M = 4) moment is 
derived for future reference. 



(c) Phase averages with scalar-valued variables of integration 

The physical moments in Table [1] can be related to phase averages. The five 
integrals with the scalar-valued variable of integration d£ will be considered first. 
With a change of variables, the moment of inertia T x (t) defined in Table [T] becomes 



Zx(t) 



dM 



[l||x(^)|| 2 -x(^)x T (i,^)] ^dc. 



(4.16) 



As all particles have the same mass by assumption, the differential mass per differ- 
ential particle dM/dip is a constant. Since we have set the ring to be unit mass, 



M = <j> p(x)d£ = 



n dM d£ 



dip = 



77 dM , „ dM , lA 
— dip = 2n— = l (4.17) 

-7T A <P dl P 



and so dM/dip = l/(2n). Consequently (|4.16p becomes a phase average, leading to 

(4.18) 



X x (t) = I||x(t)P - x (t) x r (t) 
1 



= il||x + (0|| 2 -iRe{x + (Ox?(t)} = iRe{[crx + (0] H crx + (t)} (4.19) 

with the second line following from applying the phase-averaging result f|4. llf) sepa- 
rately to each component of the matrix, and using (|4.2I) to obtain the last equality. 
Substituting for x+(f) from (|2.4I) into (|4. 19[) . we find the moment of inertia tensor 
obtains the familiar form 



Xx(«) - \ Qx(i) 



b 2 (t) 
a 2 (t) 
a 2 {t) + b 2 (t) 



Qx(*)- 



(4.20) 



It is worth pointing out that the assumption of constant mass density as a function 
of particle label, i.e. constant dM /dip, implies that the linear mass density p(x) = 
dM /d£ of the canonical ellipse is not constant; the latter changes as particles speed 
up or slow down in their flow along the ellipse periphery. 

In the same way, integrating the angular momentum over the ellipse as defined 
in Table [1] leads to 



1 



m x (i) = x(i) x x'(t) v - - Re (x+(t) x x' + *(t)} 



(4.21) 
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after reducing the integral to a phase average, then applying (|4.11[) separately to 
each component of the vector to obtain the final expression. For the scalar- valued 
zeroth moment 4(0 we may write 



4(f) 



ixf P (x)d^=iix(f)ir 



x(f)-x(f)' 



(4.22) 



on account of the fact that the phase average of x(f) vanishes. In this form we see 
that 4(0 i s the variance of the position vector x(f) around the ellipse periphery. 
As 4(f) andX u (f) are identical in form to 4(0 andX x (f) respectively, this verifies 
five of the nine correspondences in the spectral matrix expansions (|3.8[) - (|3.10j) . 



(d) Phase averages with vector-valued variables of integration 



To verify the directed line integrals in (|3.8[l - (|3.10[l . we first make an observation 
regarding the interpretation of the Hilbert transform. The obvious interpretation 
of the Hilbert transform is that it provides the particle position for a partner or 
"ghost" particle that lags ninety degrees behind the observed particle, since we have 



Hx(t) = Re {-ix+(f)} = Re { e 



i0(t)-i7r/2 



Qx(f) 



a(t) 
-16(f) 




(4.23) 



This means that the Hilbert transform contains (or more accurately, assigns) non- 
local information about what is happening at another point in the ellipse. However, 
taking the partial derivative of (|2.1j) with respect to phase tp, one finds 



d_ 

dip 



x(t,p) = Qx(t) 



— a(t) sin((/)(f) + ip) 
6(f) cos(</>(f ) + ip) 



-Kx(f , tp) 



(4.24) 



which shows that the rate of change of the position vector with variations in phase 
is the same as negative of the Hilbert transform. (In the expression Hx(t, (p), it is 
understood that T-L acts on the time variable.) This means that the Hilbert transform 
also contains (or assigns) local information, expressing the variation in the phase- 
shifted signal vector x(f , tp) as we move from one particle to the next. 

This observation allows us to simplify the computation of directed line integrals. 
For example, the circulation of the canonical ellipse becomes 

T 

(4.25) 



T x (f) EE I U T dx = f 

JC J -IT 



after a change of variables. Substituting for -g- x(f , tp) from (I4.24j) then leads to a 
phase average, and we obtain 

T x (f) = -2tt [x'(f)] T ^x(f) V = 7rlm{x^(f)x' + (f)} =7rtrSi(f) 



(4.26) 



after making use of (I4.12p . This verifies the trace of Si(f) in (|3.9[) . Similarly the 
integral for the normal vector n x (f) becomes, again using (|4.24[) and (|4. 12[) . 



n x (f) 



1 

2^ 



x x dx = — 

c 2tt 



d 

x(f,¥>) x —x.(t,(p)d<p 



= - x(f) x Hx(t) = Hx(t) x x(f) = - Im |x+(f) x x^_ 



(*)} (4-27) 
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which verifies the imaginary part of So(t) in (13.81) . In the same way we find 

n u (t) = ftx'(i) x x'(t) lf = - lm{x' + (t) x x£(t)} (4.28) 
,7 x (t) = heim[cix , (t)] T ci'Hx(t) lp = -herm Im j[crx + (t)] H crx+(i)| (4.29) 

thus verifying the four correspondences in (|3.8[) - (|3.10[) involving directed line inte- 
grals. This completes the proof of the unity of moments. 



5. Interpretation of instantaneous moments 

The recognition that the instantaneous spectral moments are identical to physical 
moments of an elliptical ring of particles calls for a close examination of key physical 
quantities — in particular, the circulation, angular momentum, and kinetic energy — 
and their relationships to instantaneous amplitude, frequency, and bandwidth. 



(a) Instantaneous frequency and circulation 

The most striking result from the unity of moments concerns the trace of the 
first-order matrix Si(t). Comparing p.5[) with (|3.9p we find 

i-T x (i) = 4(t)« x (<) = ilm{xf (t)x' + (t)} (5.1) 

which shows that the circulation of the canonical ellipse is identical to the joint 
instantaneous frequency weighted by the squared signal amplitude. This is signif- 
icant because it unifies two key quantities, the instantaneous frequency and the 
circulation, that are in use by two largely disparate communities. On the one hand, 
it provides new insight into what the instantaneous frequency is; on the other, it in- 
dicates how the time- varying circulation may be directly estimated by the analytic 
signal method. It follows that the global mean frequency is given by 

1 f°° 1 r°° 1 

uJ x = —r wtrS x Hdw = — / -r x (i)dt (5.2) 

27TC X Jo C-x J-oo n 

and therefore is interpreted as 1/n times the time-integrated circulation of the 
canonical ellipse. 

It is relevant to mention that for the canonical ellipse, ZU X , and thus the time- 
averaged circulation, must always be positive. While circulation can be either posi- 
tive or negative in general, here we have made an implicit coordinate system choice 
regarding the sign of the vertical axis with respect to the plane containing the el- 
lipse. A negative circulation can be transformed into a positive circulation through 
a three-dimensional rotation that changes the sign of the vertical axis. It is possible, 
however, that w x (t) or r x (£) may locally reverse signs and become negative. 

If we imagine that there is a two-dimensional flow in the entire area bounded 
by the canonical ellipse, then by Stokes' theorem, the canonical circulation T x (t) is 
equivalent to the vorticity integrated over this elliptical area. The spatially-averaged 
vorticity Cx(i) is given by 

C ( t )- Fx(t) -u m ° 2 (*) + b2 (*) _ 2 ^ft) f53l 
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where A(t) = ° 2 |*j + ^ 2 |*j is a measure of the ellipse shape, which, like the eccentricity, 
varies between zero for circular motion and unity for linear motion. Observe that 
for a circular vortex, X(t) = 0, and (I5.3[) reduces to Cx(i) = 2w x (<). This echoes an 
elementary result that the vorticity of a vortex in solid-body rotation is twice its 
angular frequency; thus (|5.3[) may be seen as a generalization of this result to a link 
between the time- varying canonical ellipse and its instantaneous frequency. 

A form for the joint instantane ous frequenc y w x (i) in terms of ra tes of change of 
the ellip se parameters was found bv lLillvl (|201ll ). building on work bv lLillv and Olhede 
(2010bJ) for the bivariate case. In what follows, rates of change of the ellipse angles 
will be written uj${t) = <fi'(t), LJe(t) = 0'(t), u> a (t) = a'(t), and wg(i) = (3'(t) in or- 
der to emphasize their interpretation as frequencies. The instantaneous frequency 
cj x (£) then consists of two terms 

w x (t) = w*(t) + y/l-X 2 {t) [w 9 (t)+a; (t)cosj8(t)] (5.4) 

the first being the rate of change of phase as the particle orbits the ellipse, and the 
second a generalized precession which includes the rotation of the ellipse within its 
plane, as well as azimuthal motion of the plane itself. If we associate u! a (t) with the 
vertical direction, then ui a {t) cos j3(t) is the component of the azimuthal motion or 
"spinning" of the plane containing the ellipse that is perpendicular to that plane. 
Inserting this into (|5.ip we obtain 

1 a 2 (t) + b 2 (t) 

tt r x (t) = U ^ 1 ' w*(t) + a(t)b(t) [u 9 (t) + u a (t) cos 0(t)] (5-5) 
Ztt z 

for the circulation of the canonical ellipse in terms of the rates of change of the 
ellipse parameters. It is important to point out that while this has been derived 
within the context of the canonical ellipse, its form remains valid for any set of 
ellipse parameters, and thus (|5.5[) is the circulation for any flow of the form (|2.1I) . 



(6) The normal vector and angular momentum 

Next we examine the relationship between the angular momentum vector m x (t) 
and the normal vector n x (i) to the plane instantaneously contai ning the can onical 
ellipse. The latter is given in terms of the ellipse parameters as (|Lillvl . 120111 ) 



n x (t) = ilm{x+(t) xx;(t)} 



a(t)&(t)J 3 (a(f))JiG8(t))k 



(5.6) 



where k is the vertical unit vector; consequently 7r||n x (f)|| gives the ellipse area. At 
each moment, m x (i) is split into a part that is parallel to the normal vector n x (i), 
and a part that is perpendicular to this vector, 

m 3t (*) = ni||(i)+mj.(t)= [n£(t)m(t)] n x (f) + [m(t) - m||(t)] (5.7) 

where n x (t) is a unit-length version of the normal vector, n x (<) = n x (i)/||n x (t)||. 
With straightforward but tedious algebraic manipulations, simple expressions for 
the angular momentum in terms of the ellipse parameters may be found, as is 
accomplished in |Appcndix B The parallel component is given by 

mi|(t) 



-(t)\ 



ug(t) + uj a {t) cos j8(t) + Vl-A 2 (i)^(i) 



(5.8) 
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while the perpendicular component is found to be 



lj ^ 5 = ij,(a(t))J lW t))J,(,r/2) 



I + A(f) 



cos (20(f)) sin (20(f)) 0' 
sin (261(f)) -cos (20(f)) 




oj a (t) sin /3(f) 




(5.9) 



normalizing both quantities by the magnitude of the analytic vector ||x + (f)|| 2 for 
convenience. These expressions give the exact instantaneous angular momentum of 
a ring of particles evolving according to the general kinematic model (|3.6I) . 

We see that the two vectors m x (f) and n x (f) specify the plane instantaneously 
containing the motion in two different ways, the former involving a time derivative 
of the signal, and the latter the signal's Hilbert transform. The angular momentum 
m x (f) and the normal vector n x (f) are parallel — and thus agree about the plane 
containing the ellipse — when the ellipse lies in a fixed plane, but not in general. An 
ellipse in a plane that is variable, as specified by a changing normal vector, will 
generally have a component of the angular momentum vector lying instantaneously 
within that plane. As seen in (|5.9[) . the component of the angular momentum lying 
within the plane containing the ellipse is controlled by the frequencies w a (f) and 
uip(t) that indicate motion of the plane containing the ellipse. 

The normal vector is connected in another way to the angular momentum vector 
and also to the tensor- valued first-order moment Sx{t)- One may readily show 



m x (f)+ii<(f) = ix+(f)xx' + *(f) 
(f ) + \ cr K (f ) = \ Im { [cr x+ (f )] H cr x' + (f ) } 



(5.10) 
(5.11) 



where in the latter, the right-hand-side is recognized as an expansion of the left- 
hand-side into a Hermitian and a skew-Hermitian portion. Thus the phase-averaged 
angular momentum m x (f ) and the rate of change of normal vector n x (f) are related 
as real and imaginary parts of the same complex-valued cross-product. Similarly 
the matrix J7x(f), and the cross-product matrix associated with the rate of change 
of the normal vector, are related as the Hermitian and skew-Hermitian portions of 
the same matrix. In both cases, the term involving n x (f) does not appear in the 
instantaneous moment Si(f) and thus does not contribute to the spectral matrix. 

The parallel part of phase-averaged angular momentum vector m x (f) also ap- 
pears to be closely related to the instantaneous frequency. The parallel component 
of the angular momentum (|5.8p is very similar to the form the trivariate instanta- 
neous frequency, (|5.4[) , except that the coefficients of the orbital frequency u^,(t) and 
the generalized precession ujg(t) + w Q (f) cos /3(f) have been swapped. This exchange 
between the phase "orientation" specified by </>(f ) , and the physical orientation, has 
been discussed in a very different context. In examining the general solutions to 



a freely-evolving two-dimensional shallow water elliptical vortex flow, iHohnl (|1991 



p 397) notes that reversing the roles of the orientation angle and phase angle also 
exchanges the vortex circulation and its momentum, an effect he refers to as the 
"Dedekind duality" . This comment in fact proved to be an essential clue pointing 
towards the unity of instantaneous frequency and circulation as expressed by (|5.1[) . 
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(c) Bandwidth, angular momentum, and kinetic energy 

Finally the relationship between bandwidth, angular momentum, and kinetic 
energy is examined. We will shortly show that these three quantities are related by 



lllx'WI 



IK*) I 



r x (i) 



27T||x(t)|r 



+ 2^xW 



(i) 



(iii) 



(") 



1 




2 ' 

1 


<(*) 


2 


2vr||x(t)|r 


+ 2 


K x (t) 



x(t) x x'(i) -x(t) x x'(i) 



K*)l 



(5.12) 



and thus the phase-averaged kinetic energy of the canonical ellipse can be parti- 
tioned into three distinct portions: (i) the kinetic energy associated with the circu- 
lation; (ii) kinetic energy due to expansion and contraction of the ellipse; and (iii) 
kinetic energy associated with deviation of the angular momentum vector from its 
mean value around the ellipse periphery — the angular momentum variance. This 
result remains valid for any flow of the general kinematic form (|2.ip . For the canon- 
ical ellipse, it is also seen that one-half the squared instantaneous bandwidth can be 
interp reted as an excess kinetic energy above that associated with the circulation. 

In iLillvl ([201 II ). the squared instantaneous bandwidth is found to contain four 
non-negative terms involving various rates of change of the ellipse geometry, 



<(*) 



K x (i) 
1 



|A'(0| 2 
41-A 2 (i) 



\ 2 {t) [uj e {t) +u a (t) cos f3(t)} 



(5.13) 



From left to right, these represent the effects of the rate of amplitude modulation 
K x (t), the deformation rate A'(t), and the generalized precession, together with a 
fourth term that is associated entirely with three-dimensional effects due to the 
time variation of the orientation of the plane containing the ellipse. Comparison of 
(15.131) with (|5.12[) shows that all terms except the amplitude modulation term are 
associated with the angular momentum variance. 

To prove (|5.12[) . we first note that the variance of a cross product along the 
canonical ellipse is given by, as shown in |Appendix A 



Re f (£) x Re g(t) - Re f (t) x Re g(t) v = - \\i{t) 

8 



(5.14) 



where f(t) and g(i) are analytic vectors equal to some constant times x + (t) or one 
of its time derivatives; note carefully that g(t) on the RHS is not conjugated. The 
phase variance of the angular momentum can therefore be expressed as 



x(i) x x'(t) - x(i) x x'(t) 



-(i)xx',(t)| 



(5.15) 



Article submitted to Royal Society 



22 



J. M. Lilly 



in terms of the analytic signal x+(t) and its first derivative. Lagrange's identity for 
complex- valued vectors is, as discussed in |Appendix A| 

||f xg|| 2 = ||f|| 2 ||g|| 2 -|g ff f| 2 (5.16) 

and consequently the cross-product term in (|5.15l) becomes 



|x + (t)xx' + (t)|| 2 ||x' + (i)|| 2 



]x + (t)|| 4 ||x + (*)|| 2 



<(t) 



«x(*) 



2 



«x(*) 



(5.17) 



after making use of the definition of the joint instantaneous bandwidth (|2.1ip . 
Combining (|5T7|) with (|5~T5]) then gives ([57l2| . 

6. Conclusions 

The purpose of this paper has been to introduce and interpret matrix-valued gen- 
eralizations, appropriate to the analysis of modulated oscillations in two or three 
dimensions, of the signal processing concepts of instantaneous amplitude, frequency, 
and bandwidth. These matrix-valued instantaneous moments are found to be iden- 
tical to the physical moments of the canonical ellipse that is assigned to the signal 
by the action of taking its analytic part. This gives an illuminating new interpreta- 
tion of the spectral matrix as being composed of time- varying contributions from 
the physical moments of the canonical ellipse. Since the spectral matrix is funda- 
mental to data analysis in many fields, these ideas to informatively decompose its 
variability across time could potentially find widespread application. A new result 
partitioning the kinetic energy of an elliptical vortex into three portions — associated 
with the circulation, isotropic expansion and contraction, and the variance of the 
angular momentum around the ellipse periphery — was also presented. 

A direct connection to a recently developed method for treating modulated os- 
cillations in real-world data should be mentioned. In practice, oscillatory signals 
exist not in isolation but juxtaposed with noise or other signal variability. One 
is therefore faced with the problem of extracting the presumed oscillatory signals 
from a noisy background. A powerful solution to this probl e m is w avelet ridge anal- 
ysis (|Delprat et all 119921 iMallatl . Il999t iLillv and Olhedd . l2010ah . which is based 



upon the use of an analytic wavelet transform to isolate the oscillatory signal from 
surrounding variability, leading to estimates of the analytic signal and associated 
instant aneous moments. Recent ly this method was extended to the multivariate 
case by ILillv and Olhedel ( 20121) , who also solve for the deterministic error terms 



arising from non-negligible modulation. The results of this multivariate wavelet 
ridge analysis can be used to give estimates of the instantaneous moment matrices 
corresponding to an oscillatory signal isolated from the surrounding variability, to 
which the ideas in the present paper could then be immediately applied. 

This work was support by grant #1031002 from the Physical Oceanography program of 
the United States National Science Foundation. 

Appendix A. A quartic phase average 

In this appendix we prove a phase averaging result for the quartic (M = 4) case, 
similar to the quadratic results (14. lip and (|4.12p . Observe from (|4.14p that the 
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phase average of a product of four real parts, with £„(i) = e lc ^^[c n (t) + id n (t)], is 



Re £1 (*) Re 6 (i) Re £3 (i) Re £ 4 (*) 

= - [ Cl {t)c 2 (t) + di(i)d2(t)] [c 3 (t)c 4 (t) + d 3 (t)d 4 (t)] 

+ - [ci(t)c 2 («) - di(i)d2(t)] [c 3 {t)c 4 (t) - d 3 {t)d 4 {t)} 

+ - [ Cl (t)d 2 (t) + c 2 {t)dx{t)] \c 3 {t)d 4 {t) + c 4 (t)d 3 (t)} (Al) 

o 

since all other terms involve sinusoidal functions of <f)(t) that vanish upon the phase 
integration. Comparing the last two terms on the RHS with 

Re (*)S(*)S(t)} 

= [ci{t)c 2 {t) - di [c 3 (t)c 4 (t) - d 3 (t)d 4 (t)] + 

[ci{t)d 2 {t) + c 2 {t)d 1 {t)] [c 3 (t)d 4 {t) + c A {t)d 3 {t)] (A2) 

and comparing the first term on the RHS of (|A ip with (|4.1ip . we obtain 



Re & {t) Re 6 (0 Re 6 (t) Re £ 4 (t) 



Re 6 (0 Re 6 (0 Re 6 (t) Re £ 4 (*) + g (*)& (t)« (t)} . (A3) 

A special case occurs for £ 3 (i) = £i (t) and £4 (t) — £2 (*) , in which case (|A 3|) becomes, 
after some rearranging, 



Re 6 (t) Re & (t) - Re £1 (t) Re 6 (*) 



= ^|6W| 2 |6(*)| 2 (A 4) 



which gives a simple expressions for the phase variance of Re £1 (t) Re £2 (t) in terms 
of and Comparison with (|4.1ip shows this can then be related to the 

product of the phase averages of the squares of Re£i(i) and Re&it)- 

Next we prove Lagrange's identity, presented later in (|5.16p . for two complex- 
val ued vectors f and g. This can be proven using standard vector identities (see 
e.g. Lindell . 19881 p. 4), but here we will write out components for future reference. 
The squared norm of the cross-product of f and g is given by 

l|f X gf = \fy9z ~ fz9y\ 2 + \fx9z ~ fz9x\ 2 + \fx9y ~ fy9x? 

= \fy\ 2 \9z\ 2 + \f z \ 2 \9y\ 2 + \fx\ 2 \ 9z \ 2 + \f z \ 2 \9x\ 2 + \M 2 \9y\ 2 + \f y \ 2 \9x\ 2 

- 2 Re{f y r z g* y g z } - 2 Re{f x f* z g* x g z } - 2 Re{f x f; g * x g y } (A 5) 

while the squared magnitude of the Hermitian inner product g H { is 

|g H f I 2 = \U9l + f v 9l + hg* z \ 2 = \M 2 \gx\ 2 + \f y \ 2 \g y \ 2 + \f z \ 2 \g z \ 2 

+ 2 Re{f y f: g ;g z } + 2 Re {f x f* z g* x g z } + 2 Re {Lf-glg., } . (A 6) 

Comparison of these two equations proves Lagrange's identity (|5.16p . 
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Finally we prove (|5.14l) for the phase variance of a cross product. The real-valued 
cross product Ref(t) x Reg(i) has a z component equal to 



[Ref(t) x Reg(t)], - Re f x (t) Reg y (t) - Re f x (t) Reg y {t) 



Re /„ (t) Re g x (t) - Re /„ (t) Re fo (t) ^ (A 7) 



the square of which will contain three terms, two squares and a cross term. The 
phase averages of the first two terms are 



Ref x (t)Reg y (t)-Ref x (t)Reg y (t) 



2* 1 



\f x (t)\ 2 \g y (t)\ 2 (A8) 



Ref y (t)Reg x (t)-Ref y (t)Reg x (t) 



ol/v(t)Hftc(<)l 



(A 9) 



which have simplified on account of the phase variance relation ([A 4[l . The phase 
average of the cross term is (—2) times 



Re f x (t) Re g y (t) - Re / a (t) Re ffa (t) Re /„ (t) Re g x (t) - Re /„ (t) Re (t ) 



Re (t) Re ffy (*) Re f y (t) Re 5x (i) " - Re / x (t) Re ^ (t) " Re /„ (t) Re 5x (t) " 

= 5Bfi{/ a (t) 5| ,(t)/ B *(*) 5 ;(t)} (A 10) 



using the quartic phase averaging result (|A3|) to obtain the final line. Combining 
the previous three results we find 



[Ref(i) x Reg(t)] 2 - [Ref(i) x Reg(t)] z 

= J [l/«(*)rifl^(*)r + l/*(*)| 2 lff-(*)l 2 -2Re{/ x (t) Slf (*)/-(t) ff *(t)}] (All) 



but observe that this expression occurs as one of three contributions to ||f x g|| 2 
in (I A 5[) . Considering the similar contributions from the x and y components of 
Ref(t) x Reg(t) leads to f)5 . 14[) for the phase variance of a cross product. 



Appendix B. Angular momentum derivation 

In this appendix an expression is derived for the phase-averaged angular momentum 
in terms of the ellipse parameters. Here we introduce the 2-vector 



(t) =eW*>J(0(t))r(t) =e*W 



cos 9(t) -sin6»(t) 
sm6(t) cos 9(t) 



a(t) ' 
-ib(t\ 



(Bl) 



where 3(9) is the 2x2 counterclockwise rotation matrix and r(t) = \a(t) —ib(t)] . 
The rate of change of this vector is, letting J with no argument indicate the ninety- 
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degree rotation matrix J = J(n/2), and with X(t) 

"<(*) 



_ a 2 (t)-b 2 (t) 
— a 2 (t)+b 2 (t) ' 



xV(t) 



= i</>(*) 



J(0(t)) 



«x(*) 



+ ia^(t) + Vl-A 2 (i) 



r(t) 



ws , s - 1 A' ft) 

A(t)w fl (t) + i-- 



Jr*(i) (B2) 



as may be deri ved in a few lines of algebra. Alternatively (|B 2|) follows from (69) of 
Appendix B of lLillvl (|201lh if one notes Jr(i) = - A 2 (f)r(t) + A(i)Jr(t). Recall 
that in (|B 2) , w^(f) = <j)'(t) and Wg(t) = 9'(t) as defined in the text. 

The complex- valued cross product x+(t) x x^_'(i) can be now be conveniently 
evaluated. The three- vector x + (£) defined in (|2.4p becomes 



x + (i)=J 3 Kt))J 1 (/3(i))Hx+(i) 
in terms of x+(i), where we have introduced 



H 



1 0' 
1 




(B3) 



(B4) 



as the 3x2 matrix which maps a 2-vector onto the horizontal plane in three 
dimensions. The rate of change of x + (t) is 



x',(i) = J 3 (a(t)) JiC8(t))Hx' (t) + [J 3 (a(*)) Ji(/3(t))]' Hx+(t) 



(B5) 



where the latter is found to take the form 



[J 8 (o(i))Ji09(t))]'HX f (t) = 

[w„(t)cos0(t)] J 3 (a(t))Ji(/3(t))HJ5q.(t) - {x^(t)n^(t)} n x (t). (B6) 

Observe that the real and imaginary parts of the term involving x + (i) lie within 
the plane of the ellipse, whereas the final term is collinear with the normal vector. 

With R a real orthogonal matrix such that R T = R _1 , and having unit de- 
terminant so that R is a proper rotation matrix, the cross product transforms 
as (Rf ) x (Rg) = R (f x g) for two complex-valued vectors f and g. This may 
be proven for real-valued vectors by standard vector identities, with the complex- 
valued case following by writing out the real and imaginary parts to obtain a set of 
four real- valued cross products. Also note that for two 2-vectors a and b, we have 



(Ha) x (Hb*) = [b ff Ja] k 
as may be readily verified. We then find the following expressions 



Hr(t) x Hr*(t) - [r H '(i)Jr(i)] k = iy/l - A 2 (t) ||x+(£)|| 2 k 
Hr(t) x HJr(t) = [r T (t)r(t)] k = A(t)||x+(t)|| 2 k 



(B7) 



(B8) 
(B9) 
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for the cross-product between Hr(i) and its own conjugate or itself rotated by 
ninety degrees. Combining these with (IB 5[) and (IB 6[) . we find 



[x+(t)x 
l|x + (t)|| 2 



- i 



<(*) 



Kx(t) 



+ii-^^^+w e (t)+o; a (t)cos^(t)lfi x (t) (BIO) 

for the parallel part of the complex- valued cross-product x+(t) x x+(i), i.e. that 
part parallel to the normal vector n x (i). 

To find the perpendicular part of the cross-product, we will make use of the fact 
that Ha x k = HJa for a complex two- vector a. Then 



[x + (t) x x y(t)] 



<(tR(t) 



J 3 (a(t))Ji(/9(t))H 



IM*)II 



l|x+W|| 

gives the parallel part of the cross product. Observe that this has a magnitude 
[x+(i) X x>*(t)] 



(Bll) 



|xf(0S x W| 2 



!*+(*) I 



(B12) 



and a complex direction which lies in the plane of the ellipse, rotated ninety degrees 
from the complex direction of the analytic signal. We may rearrange (|B lip to find 



[x+(t)x <•(<)]_ 



J 3 (a(t))JiC8(i))HJ 



x+(t)x?(t) 



w Q (i) sin/3(i) 



||x+(t)|P ^ w ' ||x + (t)P 

using the definition (|B 1[) of x + (i) together with 

^ c (t) = J 3 (a(t))JiC8(i))[w a (t)siii J 9(t) -^(i) 0] T 



(B13) 



(B14) 



for the rate of change of the unit normal vector. The 2x2 matrix occurring in 
(IB 13|) has a real part 



Re{x + ft)xfft)} 
l|x+(t)|| 2 



1 



a 2 (t) + b 2 (t) 



3{0(t)) 



J T (9(t)) 



a 2 (t) ' 
b 2 {t) 

cos(26>(i)) sin(20(i)) ' 
sin(20(i)) -cos(20(*)) 



(B15) 



where I2 is the 2x2 identity matrix. We then obtain (|5.9[) for the perpendicular 
part of the angular momentum vector from (|B 131) by noting HJ = J.3(7r/2)H. 
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